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Abstract. An adsorption-desorption (or parking-lot) model can reproduce qualitatively the densification 
kinetics and other features of a weakly vibrated granular material. Here we study the the two-time corre- 
lation and response functions of the model and demonstrate that their behavior is consistent with recently 
observed memory effects in granular materials. Although the densification kinetics and hysteresis are ro- 
bust properties, we show that the aging behavior of the adsorption-desorption model is different from 
other models of granular compaction. We propose an experimental test to distinguish the possible aging 
behaviors. 

PACS. 5.70.Ln,81.05.Rm,75.10.Nr 
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1 Introduction r ~ A/g where A is the maximum acceleration and g 

_ is the gravitational constant. These experiments showed 

. 5^ . The vibratory compaction of granular materials has long t^^at the density of the beads increases monotonically and 

been of importance in technological applications, but it surprisingly slowly with the number of taps, nt, for var- 

is only recently that physicists have started to investigate io^s intensities of tapping. More specifically, the density 

the process from a fundamental perspective. In particular, approaches the steady state value as 1/ In(nt). 

Knight et al. performed experiments 0| in which a column In two other studies, Nowak et al. |^J^ investigated 

containing monodisperse spherical beads was subject to a the effect of cycling F, i.e. vibrating the column for a fixed 

long sequence of taps with an intensity characterized by number of taps with a sequence of F values. They found 
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that as this parameter increases, the density of the gran- 
ular material increases. When F is subsequently reduced 
by reversing the initial sequence, the density continues 
to increase. This second branch is reversible in that it is 
retraced if the increasing sequence of F is repeated. In 
the same experiments, Nowak et al. monitored the power 
spectrum of the density fluctuations near the steady state 
for different values of F. The spectrum is distinctly non- 
Lorentzian, with the highest and lowest characteristic fre- 
quencies being separated by a non-trivial power-law-like 
regime. 

How can these effects be understood? The energy re- 
quired to displace a granular particle through a distance 
equal to its diameter is typically much larger than ksT, 
so granular materials are essentially athermal systems and 
their properties are dominated by geometrical frustration 
effects. When the system is already quite dense, even a 
small additional increase in density requires a large scale 
cooperative rearrangement of the particles. The logarith- 
mic compaction, the irreversible/reversible cycles and the 
non-Lorentzian power spectrum of the density fluctuations 
are essentially the result of this effect. 

The experimental results have stimulated a few Monte 
Carlo studies ji] and a plethora of simple models that 
incorporate geometrical frustration or quenched disorder 
(see Ref. and references therein). The latter include the 
(off-lattice) adsorption-dcsorption or parking lot model 
(PLM), Frustrated Lattice Gas (FLG) models, including 
Tetris, and one-dimensional lattice models with short range 



response properties in the parking-lot model 

dynamical constraints All these are capable of repro- 
ducing qualitatively the experimentally observed kinetics. 

Different Frustrated Lattice Gas (FLG) models |0,|,||jlO[|l3l 
have been proposed in which particles occupy the sites of a 
square lattice tilted by 45°. Each particle has two internal 
degrees of freedom which, in the Tetris realization Q, cor- 
respond to two possible orientations at a site. Neighboring 
sites can only be occupied if both particles have a favor- 
able orientation. Vibration is introduced by allowing the 
particles to move upward with a probability Pup- The con- 
nection to the compaction experiments is provided by the 
parameter 7 = —l/ln^ipup/Pdown) which plays the same 
role as F. The main virtue of this description is that it 
accounts for the layer structure of the granular material. 
In the absence of quenched disorder, the system evolves 
to a compact steady state according to a l/ln(t)-law. To 
mimic the geometric frustration of the packing, kinetic 
constraints ||7|,p^ can also be replaced by the introduction 
of quenched disorder 

Here we focus on the adsorption-desorption model in 
which hard rods are placed on a line at randomly selected 
positions with a constant rate fc_|_. If the trial particle does 
not overlap any rod already on the line, it is accepted. In 
addition, all adsorbed rods are subject to removal (desorp- 
tion) at random with a constant rate fc_ . The properties 
of the model depend only on the ratio K — k-^^/k^, with 
large values of K corresponding to small desorption rates. 
(Pioneering work was done in the limit of if — + 00 p^ , p^ ) . 
The connection to the granular compaction experiment is 
made by regarding the particles on the line as a layer in 
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the vibrating column. The effect of a tap is to eject parti- duced by the adsorption-model. All the salient features 

cles from the layer (desorption), which is followed by the of the phenomenology of vibrated granular materials are 

replacement of the particles, generally in different posi- thus qualitatively described by the model. We also show 

tions, and possibly the incorporation of additional parti- that, because of the very slow relaxation to equilibrium 

cles. 1/K plays a role similar to that of F. In this ap- (due to cooperative rearrangements of the particles), the 

proach, mechanical stability is implicitly included by the model displays out-of-equilibrium dynamical effects such 

absence of motion of the particles when the desorption is as aging jl^, and we study these properties both in the 

switch off. This microscopic and off-lattice description of response and the correlation functions. Interestingly, the 

granular compaction is suitable for describing all previ- aging behavior of the two-time correlation function of the 

ous experiments and, as we see below, is in fair agreement density fluctuations is found to be different from that ob- 

with experimental results. served for the FLG and Tetris models P,pO|. This leads us 

rr,i 1 • . • 1 -Ti • i- r_Li 1 i- to suggest an experimental way of discriminating among 

1 he kmetics and equilibrmm properties of the adsorption- oa t- j a a 

J ,. J 11 , . , . . J iH h ^ I h c l htte ^iffcrcnt models applied to vibrated granular materi- 

desorption model have been previously investigated P,P^ 15 .IqnTT: 



in the limit of small desorption rate (that is appropriate to 
describe granular compaction), the density increases very 
slowly as l/ln(t) until it nears the steady state value. 



als. 



2 Model and Simulation 



at which point the kinetics cross over to an exponen- jr ■ i ■ -i r 7 j-u i ■ -t- ■ ■ u 

^ If time IS expressed m units of k+ , the kinetics is given by 

tial form. The characteristic relaxation time Te„ of this j„ „ 

dt ^ ' K' ^ ' 

where 'P{t), the insertion probability at time t (or den- 
sity p), is the fraction of the line that is available for the 
insertion of a new particle. The presence of a relaxation 
mechanism, implies that the system eventually reaches a 
steady state that corresponds to an equilibrium configu- 
ration of hard particles with pcq = K<PcqiPcq), where pcq 
denotes the equilibrium density. At equilibrium, the inser- 
tion probability is given exactly by 

<?eq(p) = (l-p)cxp(-p/(l-p)). (2) 



final regime behaves like K'^/ln{K)'^ ||T^,[T^. Moreover, 
the density-density correlation function exhibits two char- 
acteristic time-scales: the shortest one is proportional to 
\ii{K) and is associated with a local rearrangement of par- 
ticles, whereas the longest one, Teq, corresponds to a col- 
lective rearrangement permitting an increase of the den- 
sity. The ratio of these time scales increases rapidly as K 
becomes larger, and the intermediate timescale region is 
characterized by a power law behavior. 



Here we show that the irreversible/reversible cycles ||^,^ 
and the recently demonstrated presence of memory ef- The kinetics of the model were simulated using the event- 
fects in the compaction of glass beads can be repro- driven algorithm described previously . 
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3 Densification branches and memory effects have studied the effect of an abrupt change in the des- 

orption rate 1/K on the densification kinetics. This is il- 



FoUowing the experimental procedure of compaction cycles [[^ 

lustrated in Fig. lb in which K is switched from 500 to 

we have performed simulations in which the rate of desorp- 
tion, 1//^, is first increased at a given rate and then cycled 
down and up, the simulation being stopped after the same 
time t = 40000. (Recall that 1/K plays the same role as 
the tapping intensity and that, when expressed in units of 
fc+, time is assumed to measure the number of taps; this 
leads to the reasonable behavior that the number of des- 
orption events increases as the tapping strength increases) . 
Figure|^a displays the density as 1/K varies between 10^^ Compaction cycles and memory effects characterize the 



2000 and vice-versa. The same "anomalous" behavior as 
observed experimentally by Josserand et al. ||l8|, i.e. an 
acceleration in the densification process when 1 / K is sud- 
denly lowered and the reverse phenomenon when 1/K is 
increased, is found in the adsorption-desorption model. 

4 Response functions and aging 



and 10"'^, and is similar to the behavior observed ex- response of the system to a moderate or large change of 
perimentally and numerically in the FLG models |P,p^. the tapping strength. We consider below another facet of 
Along the first (irreversible) branch the density increases the out-of-equilibrium dynamical behavior of the system 
rapidly and then passes through a maximum as 1/K in- by examining the response to a small perturbation and 
creases. When the initial sequence oil/K \s reversed {1/K the associated two-time correlation function, 
decreasing) , the density increases monotonically. When Non-equilibrium studies have been performed on FLG 
the initial (increasing) sequence of 1/K is repeated, the models: NicodemiQ demonstrated that these models ex- 
density now decreases monotonically, nearly retracing the hibit a negative response and Nicodemi and Coniglio |^ 
second branch. The residual hysteresis observed between showed that the aging regime of this model is character- 
the second and third branches is also present in other ized by C(i, t^) ^ In(f^)/ ln(i). Barrat and Loreto||l^ per- 
models|P, and it diminishes as one considers larger val- formed extensive simulations on the Random Tetris model 
ues of K and larger time intervals between changes of and clearly showed that the aging behavior is the same, 
desorption rate. Note that the densities attained in the We calculate the response function as follows. Starting 
present adsorption-desorption model are typical of a one- from an empty line, the system evolves at a fixed desorp- 
dimensional system and that more realistic values would tion rate. At the end of a fixed waiting period, t^,, two 
be obtained by employing the two-dimensional version of clones of the system are made. The original systems con- 
the model. tinues to evolve with the same desorption rate K , and the 
To further compare the model predictions to the avail- two others with K ± 5K . The response then corresponds 
able experimental data on vibrated granular materials, we to the difference in density of the copy and the original 
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system at a later time t. In order to obtain reasonable where Lw is the Lambert-W fmiction and for long times 
statistics, it is necessary to average over many indepen- the equilibrium value is attained exponentially with a char- 
dent runs (typically 10^). Since the model at equilibrium acteristic time of order / Iti{K)^ . It is worth stressing 
is equivalent to a grand canonical ensemble, the variable that with the definition adopted above of a response to 
conjugate to the density is f3fi where /3 — l/fc^T and /i In(A'), a positive response means that an increase of K 
is the chemical potential. From Eq. (|l|) at equilibrium we (i.e., a decrease in tapping strength) increases the den- 
have that /3/i = hi{K) and the response function to a sity, whereas a negative response corresponds to a density 
change in ln{K) is therefore defined as decrease. 

^(*'*-^^9M^(l)) 

Fig ^ displays two response functions for a system 

To be more precise, this is an integrated response function 

with K = 500. For a waiting time (t^ = 2000) larger 

since the perturbation is applied over an extended time. 

than the ergodic time fisl] , the system is already at equi- 

Note that when t becomes larger than the ergodic time 

librium and the response is obviously monotonically in- 

Teq , the response function attains the positive equilibrium 

creasing and in agreement with the analytical result, Eq. 

value Req — p(l— p)^. The fluctuation-dissipation theorem 

(^. At short times, the behavior of both curves is quite 

applies when > teq and we recover the time translation 

similar, but for tw — 1000 the response then decreases 

invariance R{t,tw) ~ R{t = t — t^) in which case the 

and attains a minimum before increasing again towards 

response function associated with changes in the external 

the equilibrium value (as expected in a system with in- 
field /3/j, — hi{K) (coupled to the density fluctuations) and 

terrupted aging). Fig shows the response functions in 

the density-density correlation function are related: 

the aging regime for a larger value oi K = 5000 for three 
R{t) = (5(0) — C{t) (4) different waiting times. The shape of the curves is simi- 

lar to Fig ||a, but the amplitude of the initial increase of 

where C{t) ~< 6p{T)6p{0) >. (Note that as usual for 

R{t, t^) is smaller than for K = 500. The negative mini- 
hard objects, the temperature is irrelevant and does not 

mum, that results from a competition between local rear- 
explicitly enter in Eq. (||) because it is included in \n{K) = 

rangements (leading to short-time "anomalous" response 

Pfi.) From our previous study of the fluctuations around 

discussed above) and global rearrangements (allowing for 

equilibrium ||l7|| we derive the analytical short time form 

the densification) , is less pronounced as increases. For 

which is exact for K ^ oo: 

long times, which are not accessible in the simulations, 
R{t) — - — — j-j-;^ ( I — e^^/-^^ + 0(1/ K) the response becomes positive as it approaches the equi- 

(5) librium value. 



l + Lw{K) V T + LwiK) 
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In order to investigate the out-of-equilibrium evolution of {t — tw)/tw the data for five waiting times, ranging from 
of the density-density correlation function, we consider the 25 to 5000, collapse onto a single curve. This aging behav- 
adsorption-desorption model in the limit where the des- ior, similar to that found in a large variety of systems p9|, 
orption is infinitely small||l^. In this case (K — » oo), the is thus different from the ln(t)/ln(i^) dependence of the 
relaxation time Teq is infinite and the system always FLG models. Fig. clearly demonstrates the absence of 
evolves with out-of-equilibrium dynamics. For K large, a master curve with this scaling. Moreover, no better re- 
but finite (5000), the out of equilibrium and aging be- suit is obtained when ln(i — i^)/ ln(tu,) is used in place of 
havior is interrupted since t may become larger than Teq. ln(i)/ ln(ti„). 

We have however checked that the aging behavior when This finding suggests an experimental test for decid- 
t < Teq is similar to that of the K ^ oo case. ing which phcnomcnological model, each built on different 
The normalized two-time density-density correlation assumptions: geometrical frustration in the adsorption- 
function is defined as follows: desorption model, quenched disorder or kinetic constraints 



(ju I ) — < p[t)p{tw) > < p[t) >< p[Uu) > g^^-^^ layer structure in FLG models, is more relevant to de- 

< pitwY > - < p{tw) >2 

scribe the physics of vibrated granular materials. Since all 

where the angular brackets denote an average over inde- 

predict the logarithmic compaction, the irreversible/reversible 



_ cycles, the non-Lorentzian power spectrum of the density 

Fig. m shows C{t,t^) as a function of t for different 

fluctuations and the memory effects observed experimen- 

values of the waiting time . The system exhibits aging: 

tally, a promising way to validate or invalidate them as a 

the larger t^, the longer the memory of the initial con- 
minimal theoretical approach would be to study experi- 

figuration persists. This aging regime is characterized by 

mentally the two-time correlation function of the density 

a violation of the fiuctuation-dissipation theorem and the 

fluctuations in the course of the densification process. 

breakdown of time-translational invariance [Il9|. 

When t and are large enough (but still smaller than 

J.T. thanks the National Science Foundation (CHE-9814236) 

the equilibration time), the aging behavior is usually de- 
fer financial support. P.V. thanks A. Barrat for stimulating 

scribed by a non-trivial form involving a single scaling 

discussions. 
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Fig. 1. (a) Irreversible and reversible densification branches. 
Starting from a "loose" packed state, the process consists of a 
sequence of decreasing values of K. For each K, the duration 
is 40000. The arrows indicate the way of cycle compaction. 
The dotted line is the equilibrium curve, (b) Memory effect 
at short times after a rapid change of the desorption rate: the 
full curves correspond to a process with a constant K, whereas 
the dashed curves show the kinetics of a process where K is 
switched from 500 to 2000 and vice-versa at t = 1000 



Fig. 2. Two-time response function R{t,tw) as a function of 
time t for (a) K = 500 with = 1000, 2000 and for (b) 
K = 5000 with t-u, = 10000, 20000, 30000. The horizontal lines 
correspond to the equilibrium value, Rcq = p(l — p)^ 
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ln(t) ln((t-t^^,)/t„) 



Fig. 3. Two-time density-density correlator C{t,t-u,) as 



a function of time for various waiting times tw (t- 
50, 100, 500, 1000, 5000 left to right). 




ln(t)/ln(t„) 

Fig. 4. (a) Two-time density-density correlator C{t,tw) as a 
function of (a) (t — t^^jt^ and (b) ln(t)/ In(tvK) for waiting 
times U = 50, 100, 500, 1000, 5000). 



